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We numerically simulate vortex nucleation in a BEC subject to an effective magnetic field. The 
effective magnetic field is generated from the interplay between light with a non-trivial phase struc- 
ture and the BEC, and can be shaped and controlled by appropriate modifications to the phase and 
intensity of the light. We demonstrate that the nucleation of vortices is seeded by instabilities in 
surface excitations which are coupled to by an asymmetric trapping potential (similar to the case of 
condensates subject to mechanical rotation) and show that this picture also holds when the applied 
effective magnetic field is not homogeneous. The eventual configuration of vortices in the cloud 
depends on the geometry of the applied field. 



PACS numbers: 03.75.Ss,42.50.Gy,42.50.Fx 

I. INTRODUCTION 

A spectacular property of superfiuid systems is their 
ability to support quantized vortices. These can appear 
as flux lines in superconductors to which a sufficiently 
strong magnetic field has been applied. Alternatively, in 
the case of neutral superfluids subject to sufficiently fast 
external rotation, they exist as lines of vanishing conden- 
sate density around which the velocity field flow is quan- 
tized. Both scenarios are closely interlinked, because the 
equations describing a rotating superfiuid, when studied 
in the rotating frame, mimic those of a charged superfiuid 
(a superconductor) in a magnetic field, with the Coriolis 
force playing the role of the Lorentz force. 

The dilute gas BEC is an extremely useful tool for 
probing the underlying physics of superfiuid phenomena 
because in experiment its parameters are typically much 
easier to manipulate than for other condensed matter 
systems. Moreover, because of its diluteness, it is con- 
siderably more amenable to theoretical treatment which 
more accurately reflects the experimental reality. A good 
example is in the studies of condensates rotated by an 
anisotropic potential, where the critical velocity for vor- 
tex nucleation has been found to coincide with dynam- 
ical instabilities in the surface mode excitations (TJ |2]. 
Yet, the precise mechanism for the nucleation of vortices 
is still debated, as are questions about the importance 
of thermal effects and how, if at all, the mechanism for 
vortex nucleation due to a rotating trap is related to gen- 
eration of vortices by localized stirring [3 . At the same 
time, the prospect of simulating effects related to charged 
particles in magnetic fields, means that the rotating BEC 
continues to receive a lot of attention lU |6] . 

However, rotating a condensate only provides access 



to a limited class of problems, for which the effective 
magnetic field is spatially homogeneous in the plane per- 
pendicular to the rotation axis. In addition, controlled 
stirring of a BEC can be a rather demanding task. Re- 
cent proposals to create effective magnetic fields in a 
more direct way open the door for more wide-ranging 
studies into the interaction of degenerate quantum gases 
with effective magnetic fields. One method involves us- 
ing lasers to alter the state-dependent tunneling ampli- 
tudes of atoms in an optical lattice to simulate an effec- 
tive magnetic flux [H El [9]. Another, considered here, 
exploits the interaction of A-type three-level atoms with 
two laser beams possessing relative orbital angular mo- 
mentum in an electronically induced transparency (EIT) 
configuration, such that the vector potential shows up in 
the effective equation of motion [lOl [11] for the atoms, 
which sit in a non-degenerate eigenstate of the laser-atom 
interaction. An advantage of this method is that the 
vector potential, and consequently the effective magnetic 
field, can be shaped and controlled by appropriate mod- 
ifications of the phase and intensity of the incident light. 

In this paper, we study by numerical simulation the 
influence of both homogeneous and inhomogeneous effec- 
tive magnetic fields on the dynamics of a harmonically 
trapped Bose-Einstein Condensate, and observe vortex 
nucleation for critical parameter values. The exact dy- 
namics are specific to the geometry of the trapping po- 
tential and effective magnetic field, but the the existence 
of unstable modes in the spectrum of elementary excita- 
tions as a precursor to vortex nucleation is a universal 
feature for all cases considered. Recent advances in light 
beam shaping technology, using for instance spatial light 
modulators, mean that all the potentials we consider can 
realistically be created in the laboratory [12 . 
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FIG. 1: Method for creating effective magnetic fields in de- 
generate atomic gases. On the left is the level scheme for 
the A-type atoms interacting with the resonant probe beam 
and control beam Qc- On the right is a schematic repre- 
sentation of the experimental setup with the two light beams 
incident on the cloud of atoms. The effective magnetic field 
is generated if there is relative angular momentum between 
the beams. This will occur, for example, if the probe field is 
of the form Qp ~ e^^^ , where each probe photon carries an 
orbital angular momentum M along the propagation axis z, 
and Qc is independent of the azimuthal angle. 



The paper is organized as follows: in section [ll] we out- 
line the method for creating the effective magnetic field. 



In section III we study the dynamics when a small el- 
lipticity is introduced into the trapping potential, draw- 
ing attention to key differences between the homogeneous 
and inhomogeneous effective magnetic fields. In section 
IV| we observe vortex nucleation traps of higher than 
two-fold symmetry, where the trap symmetry couples to 
higher order surface modes. Finally, in section|V|we sum- 
marize the main results. 



II. THE MODEL 

Let us briefly describe the key ingredients for creating 
an effective magnetic field using light with orbital angu- 
lar momentum (a full account can be found in [11 ). We 
consider a typical Electromagnetically Induced Trans- 
parency (EIT) setup with A-type three level atoms, char- 
acterized by two hyperfine ground levels |1) and |2) and 
an electronically excited level |3), interacting with reso- 
nant laser beams which are both propagating in the z- 
direction (see figure 1). The atoms are assumed to be 
trapped in the x — y plane with motion in the z-direction 
effectively frozen out {uOz ^ [26 . The probe beam, 
with coupling strength Qp^ and control beam, with cou- 
pling strength Qc drive the transitions |1) |2) and 
|1) 1 3) respectively. These absorption paths inter- 
fere destructively to suppress transitions to level |3) and 

/ where 

and S is the relative phase between 

the probe and control beam. Provided that gradient of 
the phase S is non zero, then, because of the spatial 
dependence of the dark state, two additional geometric 
potentials appear in the effective equation for the wave- 
function of the dark state atoms. One of these acts as 



drive the atoms to the dark state \D) — 

c 



an effective vector potential A (or Berry connection fT3]) 
and the other is a scalar potential. 

One way to achieve a non-zero phase gradient in the 
ratio of the control and probe beam coupling strengths 
is to use light with orbital angular momentum. A 
Laguerre- Gaussian beam, for example, has transverse 
field Up£{r^(j)) oc rl^lLjf' (r^) exp(— r^) exp(z-^^), where 
is the associated Laguerre polynomial, such that each 
photon carries angular momentum ih around the z-axis 
[14]. If one or both of the control and probe beam are 
of this type we may consider situations where the dimen- 
sionless parameter ( is of the form 



1^+1 

C = V^{-^) ' exp(z^0). 



(1) 



where the parameter is the intensity ratio ai r = R. 
The mean field equation for the wavefunction of the dark 
state atoms is then 

= ^(^^V + A)2* + (u + 4>)^ + g 1*1' (2) 

where, if |CP << 1, the induced vector potential is well 
approximated by [15j 



A = — — ao 



R 



which corresponds to an effective magnetic field 
haoi{iy + 1) fry 



B 



/ r \ ^-1 



(3) 



(4) 



In the << 1 limit, the scalar potential is 



v + 1 



(5) 

with the external trapping potential U for the dark state 
atoms given by 



(6) 



where Vj is the trapping potential for the atoms in 
hyperfine state j {j = 1,2). We assume that the 
(two-dimensional) coupling constant g characterizing the 
strength of atom-atom interactions within and between 
levels |1) and |2) is constant throughout the cloud which 
is valid if the inter- and intra-species s-wave scattering 
lengths are approximately equal or if |C| << 1 [15^ . 

A consequence of choosing light beams of this form is 
that the effective vector potential points in the e^ direc- 
tion, inducing a rotational motion. To elucidate the con- 
nection with rotating condensates, we use the fact that 
V • A = to re- write Eq. ^ as 



ih 



2M 



ih 
M 



A • V U, (7) 
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where V(r) = V -\- An equation of identical form 

to Eq. ^ can be obtained by considering a BEC sub- 
ject to a time-dependent trapping potential which is ro- 
tating at angular velocity 1], and studying the resultant 
Gross-Pitaevskii equation in the rotating frame in which 
the trapping potential is independent of time. The GP 
equation in this frame is as Eq. ([7| where the effective 
vector potential is given by A = Mfl x r, which would 
correspond to an effective magnetic field which is homo- 
geneous in the z-direction B = 2MQez. 

The ability to shape the phase and intensity of the 
incident light, however, means that we are no longer re- 
stricted to considering only homogeneous effective mag- 
netic fields, unlike the case of the condensate subject to 
external rotation. Actually, we can consider Eqs. ^ and 
^ as analogous to a rotation of ^ . A shearing, or differ- 
ential rotation (one at which different parts of the cloud 
are rotating at different angular velocities) is obtained if 
the vector potential A corresponds to an inhomogenous 
magnetic field. 

The possibility of studying inhomogenous effective 
magnetic fields raises many interesting questions, includ- 
ing: do we still observe vortex nucleation above a critical 
field strength? Does the same mechanism for vortex nu- 
cleation - attributed to instabilities of the surface mode 
excitations - still hold for effective magnetic fields which 
are inhomogeneous? If so what does the lattice configu- 
ration look like? 

We shall address these questions in turn, using results 
of numerical simulations of the real-time dynamics of the 
condensate as determined by Eq. ([7|, where the vector 
potential is of the form given by Eq. We shah con- 
sider the trapping potential as it appears in Eq. ([7| to 
be of the form: 



2 ^2 



V = -Mujir 



(8) 



that is, harmonic plus a small perturbation which is 
not radially symmetric. A harmonic potential can be 
achieved by a judicious choice of the external trapping 
potential (Eq. [6| to counteract the additional scalar po- 
tentials due to the light-atom interaction. Our motiva- 
tion for considering potentials of this (predominately har- 
monic) form is two- fold. First, this is of the same form 
as the rotating frame potential of atoms in a rotating el- 
liptical harmonic trap, allowing for a clear comparison. 
Also, because we are primarily interested in vortex nu- 
cleation, it is instructive to eliminate the scalar potential 
which can preclude vortex nucleation in inhomogeneous 
magnetic fields by preventing the necessary surface mode 
instabilities from occurring [15 . 

To simulate the evolution of the condensate according 
to Eq. ([7|), we first re-scale length and time variables as 



{x,y,R} 



2/i 



{x,y,R} 



respectively, where /i is the chemical potential. This al- 
lows us to write Eq. ^ in the dimensionless form 



i-fdi^l) 



(10) 



where g 



and the parameter 7 



< 1 in the 



Thomas-Fermi regime. We obtain the initial state by 
relaxing Eq. (10) in imaginary time to find the ground 



state where we constrain il) to be real, and then propagate 



this solution according to Eq. (10) by an easily imple- 
mented leap-frog method [16], checking that the energy 
and the norm are preserved throughout. This situation 
corresponds to a sudden switch on of the effective mag- 
netic field at time t = 0. 

To simulate experimental imperfections and to break 
the unrealistic levels of symmetry [17 , we add small fluc- 
tuations (^ 0.01% (F) <C grid spacing) to the trap coor- 
dinates periodically (not more often than every 10^ time 
steps). This does not change the dynamics qualitatively 
but speeds up the symmetry breaking, which would oth- 
erwise occur due to growth of numerical noise [17|. 



III. ELLIPTIC PERTURBATION TO THE TRAP 

A. Homogenous Magnetic Field {y — X) 

An already well-studied problem in respect of its equiv- 
alence to rotating condensate experiments is the homo- 
geneous effective magnetic field (i^ = 1) applied to a con- 
densate in a harmonic trap with a small elliptic pertur- 
bation [H El [H [19] 



AK, 



(11) 



(9) 



where ^Vasym. is written in the dimensionless units (Eq. 
|9|. The addition of the asymmetry in the potential dis- 
torts the flow J = ^ (V^*VV^ - V^VV^*) - A|V^p, leading to 
shape oscillations. Because the trap has a two- fold sym- 
metry, it naturally couples to the quadrupole oscillations 
which carry m = 2 units of angular momentum. At a crit- 
ical value of ao^, the frequency of the quadrupole mode 
tends to zero, and large amplitude oscillations take place. 
The critical value of a^i corresponds to an energetic in- 
stability of the quadrupole mode, occurring when the cy- 
clotron frequency associated with the effective magnetic 
field matches the bare quadrupole frequency in absence 
of the field. The large amplitude oscillations are dynam- 
ically unstable, eventually leading to nucleation of vor- 
tices in the bulk of the condensate, as shown in figure 
2. 
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FIG. 2: Numerical simulation of the density of a BEC in a 
slightly elliptic trap for 7^ = 0.001, e = 0.05 and a^l = 23, 
where the time is given in units of The almost circu- 

lar initial state becomes very elongated shortly after applica- 
tion of the effective magnetic field with cyclotron frequency 
resonant with the frequency of the quadrupole mode. The 
dynamical instability causes fragmentation of the cloud and 
subsequent nucleation of vortices. These vortices eventually 
settle into a stable lattice which is non-rotating in the labora- 
tory frame. The nucleation dynamics are in good agreement 
with the model presented in ref [17 , which considers a sudden 
switch on of rotation. 



B. Inhomogeneous Magnetic Field {y > 1) 

It is natural to contemplate whether a similar situation 
should occur for inhomogenous effective magnetic fields, 
that is, do we observe vortex nucleation due to coupling 
of the elliptic trap to a surface mode which becomes un- 
stable. An obvious place to look is at the magnetic field 
strength a^i for which the m = 2 mode becomes un- 
stable. We have previously calculated the surface mode 
frequencies of a harmonically trapped condensate subject 
to a vector potential with > 2, by linearizing around 
the non-vortex ground state solution of a BEC in ab- 
sence of effective magnetic field [15]. For ^ = 2 in Eqs. 
(|3|-(|6| and (10), the m = 2 surface mode becomes un- 
stable at aoi :^ 26 for the parameter values 7^ = 0.001, 
R = 1 and g = 1. Simulations run with this vector poten- 
tial indicate again a strong shape deformation after the 
effective magnetic field is switched on but rather than 
fragmentation followed by vortex nucleation and even- 
tual crystallization of the vortex lattice, we observe a 
more dramatic shedding of density from the cloud. The 
reason for these explosive dynamics can be explained by 
classical arguments: the effective scalar potential 'seen' 
by the rotating condensate is of the form 



V{r) = (l + e)x' + (l-e)^' + AV; 



■(7ao^)V^ (12) 



asym 



and reaches a maximum at radius rr, 



_ l+esin(26>) 

where is the polar angle. If the chemical potential /i 
1) exceeds the maximum potential energy of the trap 
V{rjnax) many atoms are expelled from the cloud. In 
such cases the system no longer becomes manageable for 
either numerical studies or actual experiments. 

This phenomena is primarily a result of our choice of 
trapping potential (Eq. (|8|), chosen to allow compar- 
isons with previous experiments. Indeed, a similar prob- 



lem exists in rotating trap experiments when the rotation 
frequency Vt exceeds the trapping frequency Ui)±_ causing 
the cloud to explode due to centrifugal effects. This has 
been overcome by adding an additional, steeper poten- 
tial to keep the atoms confined. In our case, a sufficiently 
steep confining potential exists due to the interaction of 
the atoms with the light (Eq. |5|, but the trade-off from 
using this potential in its unaltered form is that the sur- 
face modes are energetically stable for all magnetic fields 
z/ > 2, hence no dynamic instability occurs to seed vortex 
nucleation. Instead, in the next section we consider the 
resonant excitation of surface modes higher than m = 2 
as a route to achieving vortex nucleation in inhomoge- 
neous effective magnetic fields. 



IV. TRAPS OF HIGHER SYMMETRY 

Although most studies to date have considered the ro- 
tating elliptically deformed trap as a means to generate 
vortices due to the quadrupolar instability, it has also 
been demonstrated that vortex formation occurs close to 
predicted values when using a stirring potential which 
resonantly excites modes of higher multipolarity. In the 
case of a BEC in a light induced effective magnetic field, a 
trapping potential which is of rotational symmetry j, say, 
should couple to surface modes of multipolarity m = j. 
A trap with rotational symmetry higher than 2 can be 
achieved by applying a laser with multiple beam patterns 
to create the optical potential of desired symmetry [20] . 

Let us assume the asymmetric perturbation to the ra- 
dially symmetric harmonic trap (Eq. Q), written in 
dimensionless units, is of the form 



I cos ( 



I sin ( 



(13) 



where j gives the rotational symmetry, (j) is the polar 
angle and e is a dimensionless parameter. We choose the 
parameters ao£ and j to resonantly excite a particular 
mode. In figure 3 we plot the location of the surface 
mode energetic instabilities for three different effective 
magnetic fields. These instabilities were calculated in ref 
[TF, for a BEC in a radially symmetric harmonic trap. 
The addition of the potential in Eq. (13) will shift the 



mode frequencies but not by much provided e is small. 

In our simulations, we observe vortex nucleation for 
homogeneous and inhomogenous effective magnetic fields 
close to the instability of the mode coupled to by Eq. 
I\13l The cases we consider are when ^ = § and u = 2, 
wmch both correspond to an effective magnetic field 
which increases with radius r (see Eq. Q). For fields 
of this form, it is necessary to examine vortex generation 
around instabilities for surface modes higher than m = 2 
to avoid expulsion of a significant proportion of atoms. 
One should remember though that the linearization of 
the Gross-Pitaevskii equation used to calculate the ex- 
citations is less accurate for higher m modes. In figure 



Surface mode energetic instabilities for a BEC in an effective magnetic field 
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FIG. 3: Location of the energetic instability of the surface 
modes m for three different vector potentials of the form 
A ^ r^Qcf). Above the critical field strength ao^, the mode 
excitation energy is negative, so that in the presence of dis- 
sipation the system can lower its energy by going into an 
'anomolous' mode. This instability should not be confused 
with the dynamical instability whose signature would be a 
complex excitation energy. However, energetic instability is a 
prerequisite for dynamical instability. 



4, we demonstrate nucleation of vortices for the ^ = | 
field with an asymmetric perturbation to the harmonic 
trap (Eq. ([l3|) of rotational symmetry j = 4, where 
aoi is chosen to be slightly above the critical value as 
determined by figure 3. 

The presence of energetic instability is a necessary, but 
insufficient criterion for vortex nucleation; it reflects only 
that the wave-function is no longer a local minimum of 
the energy. The nucleation of vortices is due to the dy- 
namical instability, whose signature is a complex exci- 
tation energy, allowing small amplitude perturbations to 
grow exponentially in time causing a significant change to 
the initial state. Consistent with this picture, our studies 
indicate the existence of an upper as well as lower bound 
of aoi for vortex nucleation to occur. There is no upper 
bound for the energetic instability, and so the nucleation 
of vortices must correspond to a window of dynamical in- 
stability. However, because energetic instability implies 
dynamical stability [21 for the conditions outlined in this 
paper, it is still a prerequisite for the vortex nucleation. 
Moreover, the energetic stability drastically suppresses 
the energy barrier which otherwise prevents vortices en- 
tering into the cloud [22 . 

Similar to studies of refs jl9l [22l |23], we note that 
the width of the unstable region increases with e, the 
anisotropy of the trap perturbation (Eq.([l3|). For the 
parameters u = 2 and j = 5, we observe vortex nucleation 
within in the range aoi = 15.25 — 16.25 for e = 0.05 and 
ao^ = 15.0 -17.0 for e = 0.1. 

Let us also briefly mention variation in the structure 
of the vortex lattice as a function of the geometry of 
the effective magnetic field. As is well known, the con- 
densate in a homogenous effective magnetic field with a 
large number of vortices favours an Abrikosov type lat- 
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FIG. 4: Snapshots of the Density (top) and phase (bottom) 
for an inhomogeneous effective magnetic field with zy = |, 
aoi = 18.2, j = 4, where the time t/oo± = 0,20,360,700 
increasing from left to right. The m = 4 octopole surface 
mode is resonantly excited, and vortex nucleation is enabled 
by a dynamical instability. 



tice [24]. Even with a small number of vortices, we see 
that the lattice approaches a regular homogenous struc- 
ture in figure 2 for = 1. For inhomogeneous effective 
magnetic fields with > 1, in the examples studied we 
note that the vortices tend to be concentrated toward 
the edges of the cloud where the effective magnetic field 
is stronger, an effect which is enhanced as u is increased 
(see figure 5). Some additional insight into this effect 
can be gained by using the concept of diffused vorticity 
[25] . which is applicable if the there are a large number 
of vortices in the cloud, and consequently a dense vor- 
tex lattice (much more so than in the examples shown 
here). In these conditions, the course grained average of 
the velocity field Vcg associated with the lattice will be a 
smooth function of the radius r. A dense regular array of 
equally spaced vortices would thus mimic solid body ro- 
tation. Approximating the wavefunction as = y/pe**^, 
where p and S are coarse-grained averages of the density 
and phase, respectively, such that Vcg = VS*, we obtain 






FIG. 5: Density at t=450 for BEC in an asymmetric trap 
subject to the effective magnetic field with (a) = 1 (homo- 
geneous), (b) = I and (c) zy = 2. The parameters used in 
(a)-(c) were e = 0.1, j = 7, aoi = 15.5 and 7^ = 0.01. The 
nucleation of vortices is driven by resonant excitation of the 
m = 7 mode, for which the critical aoi is just below 15 for 
each of magnetic fields considered here (see figure 3). 
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an expression for the energy. 
E = [ dr 



A 

M 



Minimizing with respect to Vcg we find 



'eg 



A 



(14) 



(15) 



independent of the local density p(r), and the average 
vorticity is given by 



B 

Vxv., = -. 



(16) 



Thus, it is energeticahy favourable for the induced vor- 
ticity to mimic the effective magnetic field, implying that 
for a field of the form B (x (z^ > 1) there should be a 
higher concentration of vortices at a larger radius, as we 
indeed see in figure 5. Of course, this analysis takes no 
account of the vortex-vortex interactions (repulsive for 
vortices of the same 'charge'), which become more im- 
portant as the vortex density is increased and will tend 
to restore regularity to the lattice. 



V. CONCLUSIONS 

We have demonstrated vortex nucleation in BECs in- 
teracting with light which has a non-trivial phase struc- 



ture, where the light-BEC interaction plays the role of 
an effective magnetic field acting on the atoms. The field 
can be shaped and controlled by appropriate modifica- 
tions to the light phase and intensity, allowing us to also 
examine the role of inhomogenous magnetic fields. 

In the examples considered, we find that the nucleation 
of vortices is crucially dependent on instabilities in the 
spectrum of elementary excitations, as is known to be the 
case for condensates subject to external rotation. The 
form of the magnetic fields considered here dictates that 
the dynamics are primarily determined by the surface ex- 
citations. It is interesting to speculate whether investiga- 
tion of vortex generation due to more localized effective 
magnetic fields could bridge a connection in the theory 
between experiments where the condensate is stirred by 
an optical 'spoon' to create vortices [3] as opposed to 
rotating bucket type configurations [2|. 

The use of light-induced effective magnetic fields pro- 
vides a useful means of rotating a condensate without 
mechanical stirring. Manipulations of the trapping po- 
tential and the incident light allow a high level of control 
over the vortex dynamics, including the number of vor- 
tices in the cloud and the structure of the lattice. 
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